Spin Hall effect in a Kagome lattice driven by Rashba spin-orbit interaction 
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Using four-terminal Landauer-Biittiker formalism and Green's function technique, in this present 
paper, we calculate numerically spin Hall conductance (SHC) and longitudinal conductance of a 
finite size Kagome lattice with Rashba spin-orbit (SO) interaction both in presence and absence of 
external magnetic flux in clean limit. In the absence of magnetic flux, we observe that depending 
on the Fermi surface topology of the system SHC changes its sign at certain values of Fermi energy. 
Unlike the infinite system (where SHC is a universal constant ±^r), here SHC depends on the 
external parameters like SO coupling strength, Fermi energy, etc. We show that in the presence of 
any arbitrary magnetic flux, periodicity of the system is lost and the features of SHC tends to get 
reduced because of elastic scattering. But again at some typical values of flux {<j> = h, j, § • • •, etc.) 
the system retains its periodicity depending on its size and the features of spin Hall effect (SHE) 
reappears. Our predicted results may be useful in providing a deeper insight into the experimental 
realization of SHE in such geometries. 

PACS numbers: 73.23.-b, 72.25.Dc, 71.70.Ej 



I. INTRODUCTION 

Rapid progress in spin based information processing 
and storage devicing technologies has been metamor- 
phosed into an emerging field called 'spintronics'—, revo- 
lutionizing nanotcchnology with the plethora of concepts 
which arc particularly aimed at the exquisite control and 
manipulation of spin degree of freedom in semiconduc- 
tor structures i.e., using the spin as a career of classical 
or quantum information. Despite being a few decades 
old topic semiconductor spintronics has reignited inter- 
est to investigate the role of SO interaction in generating 
pure spin current which is the central theme of this new- 
born branch of condensed matter physics. Originating 
from the relativistic correction to the Schrodingcr equa- 
tion, SO interaction provides an all electrical way to gen- 
erate and manipulate spin current in a far precise way 
rather than the usual magnetic field based spin control. 
Longitudinal flow of unpolarized charge current through 
a sample with SO coupling can induce non-equilibrium 
spin accumulation at the lateral edges of the sample in 
transverse direction, and therefore, a pure spin current is 
established if connected through ideal leads in transverse 
direction. This is the basic phenomenon of Spin Hall 
Effect. The main source of SO coupling in mesoscopic 
systems comes from either magnetic impurities (extrin- 
sic type) or from structural or bulk inversion asymmetry 
in the confining potential of the system (intrinsic type) 
yielding Rashba or Dresselhaus type SO interaction^—. 

Few years back some theoretical proposals were made 
on the existence of intrinsic SHE in hole dopec&£ or 
electron doped 7 - semi-conducting systems where SO in- 
teraction strength is strong enough to split the Bloch en- 
ergy bands for up and down spin electrons. In this case 
a pure spin current is predicted to flow along Y direc- 
tion, in response to the longitudinal electric field along 
X direction through the infinite homogeneous system, es- 



sentially capturing the essence of a semi-classical effect. 
Sinova et al. predicted an universal value of Spin Hall 
conductivity (equal to ±g^) for two-dimensional elec- 
tron gas (2DEG) in the clean limit being independent 
of the SO coupling strength and electron density. The 
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FIG. 1: (Color online). Four-probe set-up for measuring 
spin Hall effect. An unpolarized charge current is allowed 
to pass through the longitudinal leads, 1 and 4, and because 
of Rashba SO interaction pure spin current (no charge cur- 
rent is associated with it, as they are voltage probes) flows 
through the transverse leads, 2 and 3. 

intrinsic SHE is much different from the extrinsic effect 
proposed by Hirsch^, where up and down spin electrons 
get deflected in opposite directions due to spin depen- 
dent scattering off impurities. However, the magnitude 
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of spin Hall current in the intrinsic case is expected to 
be several orders of magnitude larger than the extrinsic 
one. These theoretical anticipations have also been veri- 
fied experimentally. SHE in 2-dimensional (2D) hole^ or 
electron-^ gases have attracted a tremendous attention 
in research community. Furthermore, in some recent ex- 
periments^^ the existence of spin accumulation on the 
lateral edge of a two-terminal 2D hole gas or a 3D n-type 
semiconducting system has also been detected. 

Later on various theoretical investigations have been 
done to illustrate SHE in mesoscopic 2DEG and ring 
geometries focusing on different aspects of this intrigu- 
ing phenomena. In 2004 Sinova et al. studied SHE by 
measuring dc voltage drop in response to a longitudi- 
nal dc current^. In 2005, Sheng et al. further investi- 
gated SHE in a finite size square lattice^ to explore the 
non-quantized nature of spin Hall conductances which 
also depends on various physical parameters such as SO 
coupling strength, electronic Fermi energy, and disorder 
strength, etc. In the same year, Nikolic and co-workers 
observed^ quasi-periodic oscillations in spin Hall con- 
ductances due to spin sensitive quantum interference ef- 
fect in a finite width mesoscopic ring. In another work, 
again they studied SHE in 2DEG in detail^ both for or- 
dered and disordered cases. After that, in 2007, Nicolik 
et al. made a comparative study between extrinsic and 
intrinsic SHEs in disordered mesoscopic multi-terminal 
systems, again considering a square lattice topology^. 

Till date a wealth of literature has been formed study- 
ing SHE. But the focus of most of those works were on p- 
or n-doped semiconductors. Even the numerical calcula- 
tions for finite size systems were based on either square 
lattic o 17 ' 18 or ringi^ geometries. Later, discovery of a 
new class of spin Hall insulators^ by Murakami et al. 
establishes the fact that apart from the SO interaction 
strength the lattice structure itself plays a very signifi- 
cant role in determining spin Hall conductances through 
its band structure. In 2009, Liu et al. observed more 
exotic SHE in Kagome^i and Honeycomb^ lattices be- 
cause of their fascinating topology. But their analysis 
were based on Kubo formalism 7 -, which requires an infi- 
nite homogeneous system in clean limit. Hence a deeper 
insight into the experimental detection of such effect in 
this kind of geometry demands a quantitative prediction 
of spin Hall conductances in finite size systems. This is 
the main motivation behind our work. 

This kind of geometry is interesting not only from the 
theoretical point of view, but it has also profound ex- 
perimental significance^— . This type of lattice can be 
easily fabricated by modern patterning technique^, or 
observed in reconstructed semiconductor surfaces^ 7 -. The 
presence of a transverse magnetic field has a deeper im- 
pact on electronic properties in a kagome lattice which 
has also been reflected in our SHE study. In 1994, Nori 
et al. have studied quantum interference effect due to 
electronic motion on kagome lattice in presence of a per- 
pendicular magnetic field^ 8 -. Later, in 2002, they have 
investigated both analytically and numerically the mean 



field superconducting-normal phase boundaries in sev- 
eral 2D networks considering a transverse magnetic field. 
They have determined the transition temperature as a 
function of magnetic flux <f> passing through the small- 
est triangular plaquettes using the lattice path integral 
technique^. 

In the present work we explore different aspects of SHE 
in a Kagome lattice geometry with Rashba type SO inter- 
action in clean limit using Landauer-Buttiker formalism. 
To the best of our knowledge no such finite size spin Hall 
conductance calculation has been done for a kagome lat- 
tice geometry so far. 

The paper is organized as follows. After presenting a 
brief introduction and motivation in Section I, in Section 
II, we describe the model and theoretical formulation to 
obtain the SHC and longitudinal conductance. The nu- 
merical results are illustrated in Section III. Finally, in 
Section IV, we summarize our results. 

II. THEORETICAL FORMULATION 
A. Model and Hamiltonian 

We start describing our model which corroborates with 
the experimental set-up for observing the SHE, where 
a four-probe mesoscopic bridge is used for detection of 
pure spin current. Four ideal leads are attached to the 
central region (see Fig. [T]) which is a finite size kagome 
lattice with Rashba type SO interaction. An unpolarized 
charge current is allowed to pass through the longitudinal 
leads (lead-1 and lcad-4) inducing spin Hall current in the 
transverse direction (lead-2 and lead-3). 

A discrete lattice model is used to describe the finite 
size kagome lattice and also the side attached leads within 
the framework of tight-binding approximation assuming 
only nearest-neighbor coupling. The Hamiltonian repre- 
senting the entire system can be written as a sum of three 
terms, 

H = Hkag + Hi ea ds + H tun . (1) 

The first term represents the Hamiltonian for the finite 
size kagome lattice and it reads as, 

H kag = c l ec * + Y ''■■'■j c j + Y 1 ts ° c l P x dijh c 3 

i (ij) (ij) 

(2) 

where, 

4 = (4t c a); c ^(jl) ;e= (o °) and 

Here, e is the on-site potential energy and for a 
perfectly ordered system it is set equal to for all the 
atomic sites. c\ a and Ci. a correspond to the creation and 
annihilation operators, respectively, of an electron with 
spin a at the i-th site of the conductor, tij represents 
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the isotropic hopping strength between nearest-neighbor 
sites in presence of magnetic field. The effect of the 
magnetic field B (= V x A) is incorporated in the 
hopping term fcy through the Peierl's phase factor and 
it can be written as, 




where t is the hopping strength in absence of magnetic 
field. 4>q(= hc/e) is the elementary flux quantum. The 
specific choice of the vector potential A in this case and 
the exact calculation of the Peierl's phase factor are dis- 
cussed in the subsequent sub-section. The third term cor- 
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FIG. 2: (Color online). Schematic view of a Kagome lattice 
in presence of a perpendicular magnetic field where the struc- 
tural unit (dashed region) and the co-ordinate axes are shown. 
The length and width of the lattice strip are determined 
by two parameters, viz, N x and N v . Here, N x represents 
the number of structural units whereas N y is expressed as 
(rimax + 1) /2, where n max represents the total number of hor- 
izontal lines parallel to X axis. For this geometry n max = 3. 

responds to the Rashba type SO coupling in the system 
through spin dependent nearest-neighbor hopping term 
which introduces spin flipping in the system. The quan- 
tity t so estimates the strength of the Rashba SO interac- 
tion originated due to structural inversion asymmetry of 
the confining potential and different band offsets at the 
heterostructure quantum well interface. In this term, a 
is the spin angular momentum of the electron and d%j is 
the unit vector along the direction from i-th site to j-th 
site. 

The four metallic leads attached to the conductor are 
considered to be semi-infinite and ideal, i.e., without any 
disorder and SO interaction. The leads are described by 
a similar non-interacting single particle Hamiltonian as 
written below. 

Hleads = H a (3) 

a=l,2,3,4 



where, 

H a = ^e/c4c„ + ]P ticl n c n . (4) 

n (ran) 

Similarly, the conductor-to-lead coupling is described by 
the following Hamiltonian. 

(5) 

a=l, 2,3,4 

Here, 

Htun,a — ^c[CjC m -\- C^ n Ci] (6) 

In the above expression, t\ and U stand for the site energy 
and nearest-neighbor hopping between the sites of the 
leads. The coupling between the leads and the conductor 
is defined by the hopping integral t c . In Eq. [HI i and 
m belong to the boundary sites of the kagome ribbon 
and leads, respectively. The summation over a is due to 
incorporation of the four side-attached leads. 



B. Calculation of the Peierl's phase factor 

Now we proceed to evaluate the Peierl's phase factor 
in the term ty . 

We choose the vector potential A in the form, 

A =-Byx+^y = (-1, i 0) By. (7) 

This specific choice is followed from a literature^, and 
the purpose of doing that is solely the simplification of 
the factor J A.dl along a particular direction axis in 
this case). 

With this particular choice of A we determine Uj for 
three different types of hopping paths in the kagome lat- 
tice as follows. 

• Case 1: Our choice of gauge ensures that the compo- 
nent of A along C (see Fig. [2J axis is zero. Therefore, 
tij = t, for an electron moving along £ axis (+ve, — vc or 
its parallel direction). 

• Case 2: If we consider the motion along X axis, in 
general for the n-th line (Fig. [2]) we can write the hopping 
integral as, 

t[j = 1 6 ^° (hopping along +ve X axis) 

= t G (hopping along -vc X axis) (8) 

where, is the flux through a smallest triangle of the 
lattice, and n = 0, 1, 2, 3, . . . , (2N y - 1). 

• Case 3: Finally, we consider the hopping along fc-th 
site to i-th site and all its parallel directions (see Fig. [3]). 
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It can be shown by straightforward algebra that for an 
upward pointing triangle (Aijfc) the modified hopping 
strengths are given by, 

t k ^i = te~~ ^"+4 ) and i^ fc = k~ ln+l) (9) 
The value of n belongs to the base line of the triangle. 

Similarly, for a downward pointing triangle (Aik'j 1 ) 

downward pointing triangle 

k' j' 




j k 
upward pointing triangle 

FIG. 3: (Color online). Upward and downward pointing tri- 
angles labeled with proper site indices. 

the modified hopping integrals can be written as, 

U^y ^te-^r (n+ ^ and t k '^i 



t e -?5-( n -i)(10) 



Here, n' = n + 1, as for a downward pointing triangle the 
sites j' and k' docs not belong to the same value of n. 

Following the above prescription we incorporate the ef- 
fect of magnetic field quite easily in our lattice geometry. 
But there are other ways through which we can introduce 
the effect of magnetic field in lattice models. For exam- 
ple, Nori et al. have used a different approach of lattice 
path integral techniqu e 28 i 29 ' 31 " — to describe the effect 
of magnetic field in different lattice structures. Partic- 
ularly, this technique allows us to evaluate the physical 
quantities in terms of explicit functions for a continuously 
tunable flux. 



where, l\ and If are the charge current and spin current 
flowing through the lead-4 and lead-2, respectively. V m 
(m = 1, 2, 3 and 4) is the potential at the m-th lead. 

Now, following Landauer-Biittiker multi-probe formal- 
ism the charge and spin currents flowing through the lead 
m with potential V m , can be written in terms of spin re- 
solved transmission probabilities as^, 



T q — — \ (T aa 'V -T aa 'v ) 

1 m — r / f \ ± nm v m -'-ran v n) 



(13) 



I s = 



h ^ i^ Tnm 



)V m 4- {T nin T 7nn )V n 



(14) 

Considering linear transport regime, at absolute zero 
temperature the linear conductance (G pq ) is obtained us- 
ing Landauer conductance formula^, 



Gcrcr rpaa ( rp \ 

pq - ~T 1 pq K&F)- 



(15) 



Using Landauer conductance formula spin current 
through m-th lead can be re-written in terms of spin 
resolved conductances as, 



E r(r;H" + cit _ qU _ nil )v 
LV W nm 1 nm nm nmJ v fn 

n 

-h(G mn ~h G mn G rnn G rnn )Vn\ 



(16) 



Equation 1161 can be simplified in terms of two quantities 
defined as follows. 



mn 
ou t 
in n 



G n 

mn 



ntl 

^ ran 

*— * rrtn 



nit 

nti 



nil 
nil 

- rnm 



(17) 



Physically the term ^ n G^Vm is the total spin current 
flowing from the m-th lead with voltage V m to all other 
n leads, while the term J2 n ^mnYn defines the total spin 
current flowing into the m-th lead from the all other n 
leads having potential V n . 

Therefore, the spin current through lead m becomes, 



mY™, 



(18) 



C. Expressions of Longitudinal (Gl) and spin Hall 
conductances (Gsh) 

According to spin Hall phenomenology, in our model 
pure spin current is predicted to flow through the trans- 
verse leads, due to the flow of charge current through the 
longitudinal leads. Hence, the linear response longitu- 
dinal (4-probe) and Spin Hall conductances are defined 
as 



Gi 



J 4 



and 



G. 



H 



2e Vi - Vi 



(11) 



(12) 



Hence, following spin Hall phenomenology, in our set- 
up since the transverse leads are voltage probes, the net 
charge currents through lcad-2 and lead-3 are zero i.e., 
I\ = /| = 0. On the other hand, as the currents in the 
various leads depend only on voltage differences among 
them, we can set one of the voltages to zero without any 
loss of generality. Here, we set V4 = 0. So, from Eq. [TS] 
we have, 



{Gl 



out 
2 



Gout 
32 



Gl?) V 2 



G&Vi-G&V! (19) 



Hence, the expression of spin Hall conductance becomes 



GsH ~ Ye 



(/-tout 1 s~iout 1 y^iout \ ^2 /~ivn ^3 ,~<in 
V\ V\ 



(20) 
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This is the most general expression of SHC, but it can 
further be simplified if we assume that the leads are con- 
nected to a geometrically symmetric ordered bridge, so, 
jr = jr = 0.5. Now, in absence of external magnetic flux 
only the Rashba SO interaction does not break the time 
reversal symmetry (TRS). In general, for a time reversal 
invariant system G™ = G~° ~ a , which is equivalent to 
write G™ = —G"^. Therefore, obeying TRS the expres- 
sion for the spin Hall conductance can be written in terms 
of spin resolved transmission probabilities in a compact 
form as, 



G. 



H 



[T 4 ° 2 U * + 2T; 



out 
32 



3 T°2**] 



(21) 



One important point to be noted here is that unlike the 
charge current spin current is a vector quantity, which 
immediately gives rise to the three different components 
of spin Hall conductances (G* H , G v sH and G z sII ) and they 
are defined as follows. 



r~1X 


= IZJiVi 


-v 4 ) 




r<v 


= I y 2 /(Vx 


-v A ) 






= II/(Vi 


-v 4 ) 


(22) 



Using Eq. [21] all the three different components of SHC 
can be evaluated, only the choice of basis while construct- 
ing the matrices is important. For example, if we are 
working in cr z diagonal representation (i.e., the Z axis 
chosen to be the spin quantization axis), then Eq. [21] 
gives the ^-component of SHC. Similarly, other compo- 
nents can also be evaluated by a simple unitary trans- 
formation to the basis set. In the present work we are 
working only with the z-component of SHC. 

In a similar way from Eqs. 1111 and 1131 tlie expression of 
longitudinal conductance can be written as 



Gl 



[T 41 + 0.5 T 42 + 0.5 T 43 ] . 



(23) 



D. 



Evaluation of the Transmission Probability by 
Green's function technique 



To obtain the transmission probability of an electron 
through such a four-probe mesoscopic bridge system, we 
use Green's function formalism. Within the regime of 
coherent transport and in the absence of Coulomb inter- 
action this technique is well applied. 

The single particle Green's function operator repre- 
senting the entire system for an electron with energy E 
is defined as, 



G={E-H + irj) 1 



(24) 



where, r\ — > + . 

Following the matrix forms of H and G the problem of 
finding G in the full Hilbert space of H can be mapped 
exactly to a Green's function G^f corresponding to an 



effective Hamiltonian in the reduced Hilbert space of the 
conductor (i.e., the kagome lattice itself) and we have, 



0=Gf" = 



E — H, 



kag 



where, 



H 



t 



CJ,,H t 



(25) 



(26) 



These S a (a = 1,2,3 and 4) are the contact self-energies 
introduced to incorporate the effect of coupling of the 
conductor to the attached ideal leads. It is evident from 
Eq. [55] that the form of the self-energies are independent 
of the conductor itself through which spin transmission 
is studied. 

Following Lee and Fisher's expression for the probabil- 
ity of an electron to transmit from lead q with spin a' to 
lead p with spin a can be written as 3 -^, 



t£ = Tr\r%g r Tfg a ]. 



(27) 



r^'s are the coupling matrices representing the coupling 
between the kagome lattice and the leads, and they are 
mathematically defined by the relation, 



(28) 



Here, X£ and are the retarded and advanced self- 
energies associated with the fc-th lead, respectively 

It is shown in literature by Datta et a L 38 i 39 that the 
self-energy can be expressed as, 



iAS 



(29) 



The real part of self-energy describes the shift of the en- 
ergy levels and the imaginary part corresponds to the 
broadening of the levels. The finite imaginary part ap- 
pears due to incorporation of the semi-infinite leads hav- 
ing continuous energy spectrum. Therefore, the coupling 
matrices can easily be obtained from the self-energy ex- 
pression and is expressed in the form, 



-2Im(£; 



E. Evaluation of the Self-Energy 



(30) 



Finally, it remains the evaluation of the self-energies for 
the finite- width, multi-channel, square lattice leads. Now 
for the semi-infinite longitudinal leads (leads 1 and 4) as 
the translational invariance is preserved in X direction 
only, the wave function amplitude at any arbitrary site m 



of the leads can be written as, 
with energy 



4> m oc e %k * m * a sm(k y m y a), 



E = 2tL[cos(k x a) + cos(k y a)] 



(31) 
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In Eg. 1311 k x is continuous, while k y has discrete values 
given by, 



k y (i) 



(M+l)a 



(32) 



Here, i = 1,2,3... M. M is the total number of trans- 
verse channels in the leads, and in our case M = 2N y . 

The self-energy matrices are constructed in the reduced 
Hilbcrt space of the conductor itself. These matrices have 
non-zero elements only for the sites on the edge layer of 
the sample coupled to the leads and it is given by, 



J l(4) 



(to, n) 



M 



— siri(k y m y a)Y. r (k y ) sm.(k y n y a) 



(33) 

Y7(k y ) is the self-energy of each transverse channel and 
for a specific value of k y it becomes, 



W{k y ) 



f2 r 



2*| 



E - e(k y ) 



(E-e(k y )y 



(34) 

with e(k y ) = 2t^ cos(k y a), when the energy lies within 
the band, i.e., \E — e(k y )\ < 2th\ and 



+2 r 



E r (k y ) 



2*1 



E-e{k y ) T \(E-e(k y )Y-Atl 



(35) 



when the energy lies outside the band. Here the —ve 
sign appears for E > e(k y ) + 2|i/,| and +ve sign comes 
for E < e(k y ) - 2\t L \. 

The self-energy matrices for the other two leads (leads 
2 and 3) are also constructed in a similar way. The only 
difference in this case is that the translational invariance 
is preserved along Y direction, and accordingly, the finite 
values of k x are chosen. 



Hall conductance does not have a universal value (±^- 
as predicted in the case of an infinite 2DEG), rather it 
depends explicitly on the system parameters like Fermi 
energy (Ep), strength of Rashba SO interaction (t s0 ), 
system size, etc. The non-zero spin Hall conductance 
observed in Fig. |4] is a consequence of the fact that in 
presence of the Rashba SO interaction in the conductor, 
the up and down spin electrons flow in opposite trans- 
verse directions even in the absence of external magnetic 
flux, leading to a pure spin current, in response to the 
flow of charge current along the longitudinal direction. 

It has already been shown in literature that geometries 
like kagome lattice, graphene flakes, etc., exhibit some 
unique features of the spin Hall conductance due to their 
fascinating structures. Following Kubo formalism, Liu et 
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III. RESULTS AND DISCUSSION 

We start analyzing the numerical results by referring 
to the values of different parameters used for our calcu- 
lation. Throughout the presentation we set e = e/ = 0, 
and fix all the hopping integrals (t, ti and t c ) at the value 
1. The energy scale is measured in unit of t and choose 
the unit where c = h = e = 1. The Rashba coupling 
strength t so is also scaled in unit of t, and it is usually 
chosen as t so < t. The magnetic flux 4> is measured in 
unit of elementary flux quantum 4>o (=hc/e). 

A. Variation of spin Hall conductances with Fermi 
energy 

In Fig. @] we plot the z-component of spin Hall conduc- 
tance {G z sH ) as a function of Fermi energy (Ep) for two 
different values of Rashba coupling strengths and differ- 
ent system sizes. This figure demonstrates the fact that 
unlike the bulk or infinite system, in this case the spin 



FIG. 4: (Color online). Variation of G z sH with Fermi energy 
in the absence of external magnetic flux for two different sets 
of parameter values (shown inside the figure). 

al. have analyzed the variation of conserved spin Hall 
conductance with respect to the Fermi energy for an in- 
finite kagome lattice^. Considering a two-band approx- 
imation and treating the Rashba SO coupling as pertur- 
bation, they have analyzed the presence of various spin 
Hall plateaus (at ±^ and ±^r) with Berry phase inter- 
pretation. It has been shown that when the Fermi energy 
lies within the range — 2 < Ep < 0, the contribution to 
spin Hall conductivity due to the conventional part (a xy ) 
is —1 (in units of ^). In this case, by expanding the 
Hamiltonian around the T point upto first order in the 
Rashba coefficient t so they have established the similarity 
of the Rashba Hamiltonian with that of a semiconductor 
2DEG, and by straightforward analytical calculation spin 
Hall conductivity is obtained as — But with the in- 
crease in electron filling, within the range < Ep < 1, 
expansion of the Hamiltonian around the K point ('K- 
valley' Hamiltonian) exhibits a Dirac type spectrum with 
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linear dependence of energy on momentum, and accord- 
ingly, the value of conventional spin Hall conductivity 
gets the value 2. Furthermore the anti-symmetry at the 
band center is completely due to the particle-hole sym- 
metry of the employed tight-binding Hamiltonian. The 
variation in sign and magnitude of spin Hall conductivity 
in the range (—2 < Ep < 1) is associated with the change 
in Fermi surface topology surrounding the high symme- 
try r and K Brillouin Zone (BZ) points. The analysis 
is completely suitable for an infinite system. Here, in 
our work we try to investigate numerically whether the 
essence of these predicted features in a realistic, finite 
size, small scale, mesoscopic system are still present or 
not. 

In this figure (Fig. 0]) we observe that for a particu- 
lar system size with finite Rashba coupling strength, the 
spin Hall conductance is nearly an antisymmetric func- 
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FIG. 5: (Color online). Variation of four-probe longitudinal 
conductance as a function of Fermi energy in the absence of 
external magnetic flux. 

tion with respect to the band center i.e., G z sH = at 
Ep = 1, reflecting the particle- hole symmetry of the 
tight-binding Hamiltonian. It is a consequence of the 
fact that the spin current is defined as the difference be- 
tween spin resolved charge currents (J^ and I±), and the 
spin current carried by negatively charged electrons at 
Ep < I can be interpreted as the propagation of pos- 
itively charged holes with opposite spin in opposite di- 
rection. Apart from that, sign reversals take place also 
at Ep = and at Ep = 2 as already predicted for the 
infinite system. The oscillatory behavior in SHC for a 
smaller system size is entirely the finite size effect. The 
feature of sign reversals in SHC gets prominent with the 
increase of system size. Hence, our numerical results for 
finite sized systems provide all the essential features of 
SHC pattern observed in an infinite kagome lattice. 



addition of the two other semi-infinite transverse finite 
width leads allows some extra phase-breaking paths and 
thereby lifting the transmission zeros and hence broaden- 
ing the conductance peaks. The conductance spectrum 
reveals itself the energy eigenstates of the finite size sys- 
tem. The presence of Rashba coupling does not affect the 
energy spectrum in a significant way apart from shifting 
of the energy eigenvalues a little. 



C. Spin Hall conductance as a function of Rashba 
coupling strength 

The dependence of spin Hall conductance of a finite 
size kagome lattice on the Rashba coupling strength at 
the typical energy Ep = —0.8 is illustrated in Fig. [5] Un- 
like the results predicted by using linear response theory, 
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FIG. 6: (Color online). Z-component of spin Hall conductance 
as a function of Rashba coupling strength in the absence of 
magnetic flux for a typical system size and at a particular 
Fermi energy. 

in this case G Z H — > as t so — > 0, and G Z H shows non- 
negligible, non-zero values within the range 0.1 < t so < 5. 
The strength of Rashba SO interaction can be tuned ex- 
ternally by using a gate voltage. Although the physically 
accessible range of t so (in unit of t) is 0.001 to 1, but 
here we plot the spectrum for a wider range, a part of 
which may be beyond the experimental reach till date. 
For a sufficiently large SO coupling, G z sH again drops to 
0, because large SO coupling in the conductor forms a 
large potential barrier for the incident electrons, yielding 
a very small probability to transmit through the conduc- 
tor. 



D. Effect of Magnetic field 



B. Variation of longitudinal conductance with 

Fermi energy 1- Energy-flux characteristics 



In Fig. [5] we explore the variation of four-probe longi- 
tudinal conductance in presence of finite Rashba interac- 
tion strength (t so ^ 0) as a function of Fermi energy. It 
exhibits quite a similar behavior to our previous investi- 
gation^ of two-terminal longitudinal conductance. The 



Figure [7] depicts the energy- flux characteristics of a fi- 
nite size (N x = 6 and N y = 2) kagome lattice both in 
the presence and absence of Rashba type SO coupling. 
The energy eigenvalues are obtained by diagonalizing the 
Hamiltonian and the energy-flux spectrum is often called 
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the Hofstader spectrum. Form Fig. [7ta) it is clearly visi- 
ble that a highly degenerate level exists at E = —2t when 
(f> is set equal to 0, and the application of a very small 
non-zero flux starts to break the degeneracy. The pres- 
ence of even a very small magnetic flux affects the phase 
of the electronic wave function and thus tends to destroy 
the quantum interference and eventually breaks the flat 
band^i. Again at the half flux-quantum (<fi = 4>o/2), the 
flat band re-appears because of the interference effect but 
at the energy E = 2t. The energy levels are periodic in 
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FIG. 7: (Color online). Energy-flux characteristics of a fi- 
nite size kagome lattice in (a) absence and (b) presence of 
Rashba type SO coupling. The position of the highly degen- 
erate energy level at E — — 2 is slightly shifted due to Rashba 
interaction but the degeneracy is not broken. 

cf>, showing (f>Q (which is set to I in our calculation) flux- 
quantum periodicity and the energy spectrum is mirror 
symmetric about </> = 0.5. The presence of a small non- 
zero SO interaction does not affect the degeneracy much 
which is clearly noticed from Fig. [7th). 



2. Effect of magnetic flux on spin Hall conductance 

In Fig. [5] we plot the nature of SHC as a function 
of Fermi energy for different values of magnetic flux <f>, 
where (a), (b) and (c) correspond to <j> = 0, 0.32 and 
0.5, respectively. From the spectra it is observed that 
by applying an arbitrary magnetic flux (0 = 0.32) the 
features of SHE tend to get diminished, but interestingly 
at <j) — 4>q /2, all the above mentioned features re- appear 
being just the inverse of the case for <f> — 0, in accordance 
with the change of band structure with the application 
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FIG. 8: (Color online). Variation of G z sH as a function of 
Fermi energy for three different values of flux, where (a) (j> = 
(b) 4> = 0.32 and (c) <j> = 0.5. 



of magnetic flux. These features can be illustrated as 
follows. 

If we write the Schrodinger equations for three different 
sites in a unit cell (Fig. of an infinite kagome lattice^, 
we get, 

Ei/)(MA,NA,1) = ip(MA, NA, 2) + ip(MA, (N + I) A, 2) + 



"* {N+ i } H(M - I) A, (N + 1)A, 3) 



Et/}(MA,NA,2) 



E^}{MA,NA,Z) 



e *o 
ip(MA, NA, I) + 
ip((M-l)A, (JV- 1)4,1) + 

-i8ivN<!> 

e *o ip((M - I) A, NA, 3) + 

e lS ^ ± ip{MA,NA,2) 

e l8 ^ N * tp(MA, NA, 2) + 

i/j((M + l)A, NA, 2) + 

e 



~\Z ( N+ i~)vP(MA, NA, 1) + 
e = ^ ±( ' N ~ i ' > ip(MA, (N - 1)A, 1) 
where, ip{MA,NA,f) denotes the wave amplitude at a 



(36) 
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particular site (MA,NA,j) of the unit cell. Here, M 
indicates the M-th triangle along the X axis, N indi- 
cates N-th triangle along the ( axis and j = 1,2 and 3 
represents the vertex in the triangle. 

Here, for any arbitrary flux, translational invariance is 
lost along all directions and the hopping integrals become 
different. As a result the sharp features like the sign 
reversal about the symmetry points of the Fermi surface 
get reduced because of the elastic scattering, which is 
visible in Fig. [S^b). Very interestingly, again at <f> = 
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FIG. 9: (Color online). Unit cell configuration of an infinite 
kagome lattice. 

c/>o/2,it can be shown from the above equations that the 
translational invariance is again restored along both the 
X and £ directions which can be deducted from the above 
equations, and therefore, scattering effect gets completely 
suppressed and SHC re-appears but with inverse manner 
due to the change in band structure. 

It can also be observed that for 6 = (n and m being 

^ 8m v ° 

integers), translational invariance is retained (increasing 
the unit cell dimension) and hence system size has an 
important role as we are dealing with a finite size kagome 
lattice. In this case, some anisotropy is introduced in the 
system through the hopping terms. Therefore, the SHC 
pattern changes accordingly due to the changes in the 
Fermi surface topology. In our case (N x = N y = 8) we 
observe that for </> = j and | a regular pattern in SHC is 
obtained along with the sign reversal at the band-center 
(see Fig. [T0|) . The application of an external magnetic 
flux breaks the time reversal symmetry (TRS), but still 
SHE is still observed though reduced due to the clastic 
scattering. 

The reason is that in the presence of an external mag- 
netic field along the Z direction, an ordinary Hall volt- 
age appears along the transverse edges of the sample and 
since the transverse leads are voltage probes no charge 
current will flow through them. Therefore the tendency 
of accumulation of charges of both spins somehow tries 
to reduce the features of SHE, which is entirely an effect 
due to the spin accumulation along the transverse edges 
and the flow of spin current in that direction, if leads are 
connected. 

Before we end this section, we would like to point out 
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FIG. 10: (Color online). Variation of G z sH as a function of 
Fermi energy for (a) <j) = \ (b) <j) = |. 

that in a kagome lattice geometry Nori et al£L have de- 
scribed several physical phenomena in the presence of a 
transverse magnetic field. They have analyzed the for- 
mation of cusps in the variation of superconducting tran- 
sition temperature as a function of magnetic field in this 
geometry at </> = g, |, |, |, |, etc., considering lower or- 
der approximation in the lattice path integral method 
which directly reflects the lattice topology on the elec- 
tronic properties. Though our present analysis is quite 
different from their analysis, but here we also establish 
that the SHE exactly reappears at these typical values of 
magnetic flux <p which reflects the nature of the lattice 
model. 



IV. CLOSING REMARKS 

In conclusion, in the present paper we have studied dif- 
ferent aspects of mcsoscopic spin Hall effect induced by 
Rashba type SO interaction in a kagome lattice geometry 
attached to four finite width probes both in the presence 
and absence of external magnetic flux. We have evalu- 
ated spin Hall conductance (SHC) and longitudinal con- 
ductance for a finite size system in the clean limit using 
four-terminal Landauer-Biittikcr formalism and Green's 
function technique. In the absence of magnetic flux, we 
have observed that due to the change in Fermi surface 
topology SHC changes its sign at certain values of Fermi 
energy, along with the band center. Unlike the infinite 
system (where SHC is a universal constant ±^r), here 
SHC depends on the external parameters like SO cou- 
pling strength, Fermi energy, etc. We have shown that 
in presence of arbitrary magnetic flux, periodicity of the 
system is lost and the features of SHC get suppressed 
because of the weak elastic scattering, but not lost com- 
pletely. On the other hand, at some typical values of flux 
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(4> = i. j, I . . ., etc.) the system retains its periodicity 
again depending on the system size and the features of 
spin Hall effect (SHE) re-appears. 

It is also important to note that in our theoreti- 
cal model we have included the effect of magnetic field 



through the phase factor in the hopping term and ignored 
the Zeeman term in the Hamiltonian since in this geome- 
try, the Zeeman coupling strength (~ 0.03 meV) is much 
smaller than the hopping integral (~ 1.5 meV) which has 
been discussed in the earlier work^i. 
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